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Abstract . This paper is concerned with the optimization and 
guidance of trajectories for coplanar, aeroassisted orbital 
transfer (AOT) from high Earth orbit (HEO) to low Earth orbit 

l 

(LEO). In particular, HEO can be a geosynchronous Earth orbit 
(GEO) . It is assumed that the initial and final orbits are 
circular, that the gravitational field is central and is governed 
by the inverse square law, and that at most three impulses are 
employed, one at HEO exit, one at atmospheric exit, and one at 
LEO entry. It is also assumed that, during the atmospheric pass, 
the trajectory is controlled via the lift coefficient. The 
presence of upper and lower bounds on the lift coefficient is 
considered . 

First, optimal trajectories are computed by minimizing the 
total velocity impulse (hence, the propellant consumption) 
required for AOT transfer. Use is made of the sequential gradient- 
restoration algorithm (SGRA) for optimal control problems. It is 
shown that the optimal trajectory includes two branches: a 
relatively short descending flight branch (branch 1) and a long 
ascending flight branch (branch 2). 

In branch 1, the path inclination ranges from a few degrees 
negative to zero and is nearly a linear function of the altitude; 
in branch 2, the path inclination ranges from zero to a fraction 

of a degree positive and is a slowly varying function of the 
altitude. Velocity depletion takes place along the entire 
atmospheric trajectory, but is concentrated mostly in the terminal 
part of branch 1 and the beginning part of branch 2. As 
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the ratio of HEO radius to LEO radius increases, the minimum 
altitude of the optimal trajectory decreases, implying that a 
deeper penetration into the atmosphere is required. 

Next, attention is focused on guidance trajectories capable 
of approximating the optimal trajectories in real time, while 
retaining the essential characteristics of simplicity, ease of 
implementation, and reliability. For the atmospheric pass, a 
feedback control scheme is employed and the lift coefficient 
is adjusted according to a two-stage gamma guidance law. For 
branch 1, the gamma guidance is a linear path inclination guidance, 
for branch 2, the gamma guidance is a constant path inclination 
guidance. The switch from branch 1 guidance to branch 2 guidance 
is governed by the requirement that a specified apogee be reached, 
following the atmospheric exit. Computer simulations show that, 
by proper selection of four guidance parameters (the entry path 
inclination Yq* the target altitude h,^ of branch 1, the switch 
velocity V , and the target path inclination y T of branch 2) , the 

s 

gamma guidance trajectory can be made close to the optimal trajectory. 

Further improvements are possible via a modified gamma 
guidance, which differs from the gamma guidance as follows: in 
the gamma guidance, the parameters V^, are preselected; in 
the modified gamma guidance, V^, y^ are adjusted in flight with a 
predictor-corrector algorithm; also, the target altitude of the 
modified gamma guidance is lower than that of the gamma guidance. 
Computer simulations show that the modified gamma guidance trajectory 
is superior to the gamma guidance trajectory in the following sense. 
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it is more stable with respect to dispersion effects 
arising from navigation errors, variations of the atmospheric 
density, and uncertainties in the aerodynamic coefficients. 

A byproduct of the studies on dispersion effects 
is the following design concept. For coplanar aeroassisted orbital 
transfer, the lift-range-to-weight ratio appears to play a more important 
role than the lift-to-drag ratio. This is because the lift-range-to- 
weight ratio controls mainly the minimum altitude (hence, the 
peak heating rate) of the guidance trajectory; on the other hand, 
the lift-to-drag ratio controls mainly the duration of the 
atmospheric pass of the guidance trajectory. 

Key Words. Optimization, guidance, gamma guidance, modified 
gamma guidance, parameter dispersion effects, flight mechanics, 
astrodynamics , aeroassisted orbital transfer, sequential gradient- 
restoration algorithm. 
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1 . Introduction 

Saving propellant weight and increasing the payload are 
among the most important problems of space transportation. Orbital 
transfer from high Earth orbit (HEO) to low Earth orbit (LEO) 
can be made more economic if the aeroassisted orbital transfer 
(AOT) mode is employed. In the AOT mode, use is made of the 
aerodynamic forces in order to achieve the proper amount of 
velocity depletion during the atmospheric pass. Here, the intent 
is to achieve a specified apogee following the atmospheric 
exit, while minimizing the overall propellant consumption and 
keeping the peak heating rate within reasonable bounds during the 
atmospheric pass. 

Aeroassisted orbital transfer is not only important for 
HEO-to-LEO transfer maneuvers, but may prove to be indispensable 
for future planetary flights. In particular, this statement refers 
to lunar return vehicles, Mars exploration vehicles, and Mars 
return vehicles. 

Over the past several years, considerable research has been 
done on two aspects of coplanar, aeroassisted orbital transfer: 
trajectory optimization (Refs. 1-7) and trajectory guidance 
(Refs. 8-12) . Concerning trajectory guidance for lift-modulated 
AOT, see the work of Mease and McCreary (Ref. 9) , Lee and Grantham 
(Ref. 10), and Miele, Wang, and Lee (Ref. 12). In particular, 
in Ref. 12, a two-stage guidance scheme, consisting of the 
combination of target altitude guidance and target path inclination 
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guidance, was developed for the atmospheric pass of an AOT 
spacecraft, akin to the target altitude guidance already developed 
'for the abort landing of an aircraft in a windshear (Ref. 13). 

This paper continues the work of Ref. 12 and develops a 
two— stage gamma guidance for the atmospheric pass of an AOT 
spacecraft, akin to the gamma guidance already developed for the 
abort landing of an aircraft in a windshear (Ref. 14). Indeed, 
there are similarities between these two situations: (i) both 

trajectories are characterized by descending flight, followed 
by near-horizontal flight, followed by ascending flight; (ii) 
both trajectories are characterized by energy dissipation? 

(iii) both trajectories are dangerous unless proper guidance and 
control are applied. However, there are also differences between 
these two situations: (i) the performance index of the AOT 

spacecraft is the characteristic velocity, while the performance 
index of the abort landing aircraft is the altitude drop; (ii) 
the energy dissipation of the AOT spacecraft is due to the 
aerodynamic forces, while the energy dissipation of the abort 
landing aircraft is due to the combination of shear and downdraft; 
(iii) reaching a specified apogee after the atmospheric pass 
is essential for the AOT spacecraft, while reaching a specified 
aititucie after the windshear encounter is not essential for the 
abort landing aircraft. The similarities suggest that the gamma 
guidance scheme developed for abort landing in a windshear can 
be adapted in concept to the atmospheric pass of an AOT spacecraft. 
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The differences suggest that major modifications are necessary. 

With the above ideas in mind, this paper is organized as 
follows. Section 2 contains the notations, and Section 3 presents 
the system description. In Section 4, we study optimal trajectories, 
namely, trajectories minimizing the propellant consumption 
required for orbital transfer. We observe that, during the 
atmospheric pass, the optimal trajectory includes a relatively 
short descending flight branch (branch 1) and a long ascending 
flight branch (branch 2) . We also observe that the path inclination 
of the optimal trajectory can be approximated by a linear function 
of the altitude in branch 1 and by a constant value in branch 2. 

For the guidance trajectory, this leads to a two-stage 
gamma guidance (Section 5) , in which the guidance law is a linear 
path inclination guidance in branch 1 and a constant path inclination 
guidance in branch 2. Next, via computer simulations, we select 
the four guidance parameters of the gamma guidance, namely, the 
entry path inclination y the target altitude h^ of branch 1, 
the switch velocity V , and the target path inclination y t of 
branch 2. Having developed the gamma guidance, we then proceed 
to developing the modified gamma guidance (Section 6), in which 
the parameters y n , h are preselected, while the parameters V , 

Y are adjusted in flight via a predictor-corrector algorithm; 

T 

in addition, the target altitude of the modified gamma guidance 
is lower than that of the gamma guidance. Finally, we test the 
modified gamma guidance with respect to dispersion 
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effects ari 
atmospheric 
(Section 7) 


sing from navigation errors, variations of the 
density, and uncertainties in the aerodynamic coefficients 
The conclusions are given in Section 8 . 
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Notations 

Throughout this paper, the following notations are employed 
= drag coefficient; 

= zero-lift drag coefficient; 

= lift coefficient; 

= drag, N; 

2 

= local acceleration of gravity, m/sec ; 

= altitude, m; 

= r - r = thickness of the atmosphere, m; 
a e 

. 2 

= heating rate, W/m ; 

= induced drag factor; 

= gain coefficient; 

= lift, N; 

= ma s s , kg; 

= radial distance from the center of the Earth, m; 

= radius of the Earth, m; 

= radius of the outer edge of the atmosphere, m; 

^ 2 

= reference surface area, m ; 

= T/t = dimensionless time; 

= running time, sec; 

= velocity, m/sec; 

= /(y/r ) = circular velocity at r = r & , m/sec; 

= path inclination, rad; 

3 . 2 

= Earth's gravitational constant, m /sec ; 

3 

= air density, kg/m ; 

= final time, sec; 
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AV = characteristic velocity, m/sec. 

Subscripts 

'0 = entry into the atmosphere; 

I = exit from the atmosphere; 

00 = exit from the initial orbit; 

II = entry into the final orbit. 

Superscripts 

. = derivative with respect to dimensionless time; 

= condition following the application of a velocity impulse 
or nominal condition. 
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3 . System Description 

Xn this section, we consider coplanar , aeroassisted orbital 
transfer from high Earth orbit to low Earth orbit. 

We employ the following assumptions: (i) the initial and final 

orbits are circular; (ii) three impulses are employed, one at 
the exit from the initial orbit, one at the exit from the 
atmosphere, and one at the entry into the final orbit; and 
(iii) the gravitational field is central and is governed by the 
inverse square law. The four key points of the maneuver are these: 
point 00, exit from the initial orbit; point 0, entry into the 
atmosphere; point l,exit from the atmosphere ; and point 11, entry 
into the final orbit; see Fig. 1. 

The maneuver starts in high Earth orbit with a tangential 
propulsive burn, having characteristic velocity AVq Q , at point 00; 
here, the spacecraft enters into an elliptical transfer orbit, 
connecting the points 00 and 0; this elliptical transfer orbit 
is such that its apogee occurs at r^. At point 0, the spacecraft 
enters into the atmosphere; after traversing the upper layers 
of the atmosphere, it exits from the atmosphere at point 1; during 
the atmospheric pass, the velocity of the spacecraft is depleted, due 
to the aerodynamic drag. At point 1, the maneuver continues 
with a tangential propulsive burn, having characteristic velocity 
AV 1< * then, the spacecraft enters into an elliptical transfer orbit 
connecting the points 1 and 11; this elliptical transfer orbit is 
such that its apogee occurs at r,.* 


The maneuver ends with a 
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tangential propulsive burn, having characteristic velocity 
AVu , at point 11; here, the spacecraft enters into the low 
,Earth orbit, in that the magnitude of AV^ is such that the 
desired circularization into LEO is achieved. 

3.1. Atmospheric Pass . For the atmospheric portion of the 
trajectory of the AOT vehicle, we employ the following hypotheses: 

(i) the atmospheric pass is made with engine shut-off; hence, 

in this portion of the flight, the AOT vehicle behaves as a 
particle of constant mass; (ii) Coriolis acceleration terms 
and transport acceleration terms are neglected; (iii) the 
spacecraft is controlled via the lift coefficient; 

(iv) the aerodynamic forces are evaluated using the inertial 
velocity, rather than the relative velocity; (v) under extreme 

hypersonic conditions, the dependence of the aerodynamic coefficients 
on the Mach number and the Reynolds number is disregarded. 

3.2. Differential System. With the above assumptions, and 


upon normalizing the flight time to unity. 

the equations of motion 

are given by 


• 

h = t [V siny ] , 

(la) 

V = t [-D/ m - gsiny] , 

(lb) 

y = T [L/mV + {V/r - g/V)cosy], 

(lc) 


with 0 < t <_ 1. In the above equations, 

2 2 

= r + h, g = y/r = y/(r + h) , 

0 


r 


( 2 ) 
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where y denotes the Earth's gravitational constant. In addition, 

the aerodynamic forces are given by 

D = (1/2)C d oSV 2 , L = (l/2)C L pSV 2 , (3, 

with p = p(h). In particular, if a parabolic polar is postulated, 
the relation between the drag coefficient and the lift coefficient 
is given by 


C = C + KC ‘ 
D DO L 


(3b) 


3 . 3 . Control Constraint . To obtain realistic solutions, 
the presence of upper and lower bounds on the lift coefficient 
is necessary. Therefore, the two-sided inequality constraint 


C T < C T < C_. 
La — L — Lb 


(4) 


must be satisfied everywhere along the interval of integration. 

3.4. Boundary Conditions. At the entry into the atmosphere 
(t = 0 ) and at the exit from the atmosphere (t = 1 ) , certain 
static and dynamic boundary conditions must be satisfied. Specifically, 
at atmospheric entry, we have 


h 0 = H ' 


(5a) 


2 „ 2 2 


00 


< 2V a - V 0 } ' 2 r 00 r a V a + r a V 0 COS ^0 = °' 


where H is the thickness of atmosphere and V_^ is the circular 
velocity at r = r . In addition, at atmospheric exit, we have 

3 . 


(5b) 
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h ± = H, 

2 2 ~ 2 2 2 ~ 2 2 n 
r ll (2V a - V 1 > - 2r ll r a V a + r a V 1 cos Y 1 = °' 


(6a) 

(6b) 


Note the following: (i) the entry condition (5a) implies that 

r = r ; (ii) the entry condition (5b) arises from energy 
0 a 

conservation and angular momentum conservation applied to the 
HEO-to-entry transfer orbit; (iii) the exit condition (6a) implies 
that r = r ; (iv) the exit condition (6b) arises from energy 

1 H 

conservation and angular momentum conservation applied to the 
exit-to-LEO transfer orbit; (v) because the velocity impulse 

is applied at atmospheric exit, the following relation holds: 



(7a) 


here, V x and V 1 denote the values of the exit velocity before 
and after the application of the propellant burn; (vi) in the 

light of the exit condition (6b) , the value of the exit velocity 
a fter the application of the propellant burn can be written as 


V 1 - v a / ‘ [2(r Ir 


r ll r a> /<r H 


2 2 

- r cos Yt ) 1 

cl 


3.5. Summary . The relations governing the 
atmospheric pass include: the differential system (l)-(3), the 
control constraint (4), and the boundary conditions (5) -(6). In this 
formulation, the independent variable is the time t, 0 £ t <_ 1. 

The dependent variables include the state variables h(t) , V(t) , 
y(t), the control variable C L (t) , and the parameter x. 
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3.6. Experimental Data . The following data are used in 
the numerical experiments for optimal trajectories and guidance 
trajectories (Sections 4-7) . 

Spacecraft . For the spacecraft, it is assumed that the mass 

2 

per unit reference surface area is m/S = 300 kg/m ; the zero-lift 
drag coefficient is C dQ = 0.1; the induced drag factor is K = 1.11; 
the lift coefficient for maximum lift-to-drag ratio is C LE = 0.3; 
the maximum lift-to-drag ratio is E = 1.5; the bounds on 
the lift coefficient are = -0.9 and = +0.9. 

Physical Constants . The major physical constants used in the 
computations are as follows: the radius of the Earth is r g = 6378 km; 
the radius of the outer edge of the atmosphere is r = 6498 km; 

a. 

the thickness of the atmosphere is H = 120 km; and the Earth's 

3 2 

gravitational constant is y = 0.3986E+06 km /sec . 

Transfer Maneuvers . Three transfer maneuvers are considered, 
involving different values of the HEO radius, but the same value 
of the LEO radius. To describe these maneuvers, let a and 6 
denote the dimensionless ratios 


° = r 00 /r a ' 


6 = r ll /r a 


Case 1. The HEOl radius is r Q0 = 12996 km, a = 2. The 
LEO radius is r^ = 6558 km, 8 = 1.00923. 

Case 2. The HE02 radius is r^ = 25992 km, a = 4. The 
LEO radius is r. 1 = 6558 km, 8 = 1.00923. 

Case 3. The HE03 radius is r^^ = 42164 km, a = 6.48877. 
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The LEO radius is r^ ~ 6558 km, 6 = 1.00923. Note that HE03 = GEO. 

Atmospheric Model . The atmospheric model used is 
•the US Standard Atmosphere, 1976 (Ref. 15). In this model, the 
values of the density are tabulated at discrete altitudes. For 
intermediate altitudes, the density is computed by assuming an 
exponential fit for the function p (h) . This is equivalent to 
assuming that the atmosphere behaves isothermally between any 
two contiguous altitudes tabulated in Ref. 15. 

Heating Rate . The heating rate HR is computed with 
the relation 

HR= c/lp/p R ) (V/V R ) 3 ' 08 . (9) 

3 

Here, p = 0.39957E-02 kg/m is a reference density (density at 
R 

the reference altitude h R = 40 km) and V R = V Q = 7.832 km/sec 
is a reference velocity. The constant C represents the heating 

. 2 

rate at V = V_ and h = h D ; its value is assumed to be C = 348.7 W/cm . 
R R 
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4 . Optimal Trajectories 

4.1. Performance Index . Subject to the previous constraints, 
different AOT optimization problems can be formulated, depending 
on the performance index chosen. The resulting optimal control 
problems are either of the Bolza type or of the Chebyshev type. 

In this paper, only one performance index is considered, the 
minimum energy required for orbital transfer. A measure of this 
energy is the total characteristic velocity AV, the sum of the 
characteristic velocity AV^q associated with the propulsive burn 
from the initial orbit, the characteristic velocity AV^ associated 
with the propulsive burn at the exit of the atmosphere, and 
the characteristic velocity AV^ associated with the propulsive 
burn into the final orbit. Clearly, 


I = AV = AV 


00 


AVj- + 


AV 


11 


(10a) 


with 


AV 00 = /lr a /r 00 )V a ‘ <r a /r 00> V 0 COSY 0 ' 


Av i - v i - V 


AV 11 = /lr a /r ll )V a ‘ <r a /r ll ,V l COSY l' 


(10b) 

(10c) 

(10d) 


In the last two equations, V^ is supplied by Eq. (7b). 

4.2. Numerical Results . Optimal trajectories were computed 
by minimizing the performance index (10) , subject to the constraining 
relations. Three transfer maneuvers were considered: HEOl-to-LEO, 
HE02-to-LE0, and HE03-to-LE0; see Cases 1,2,3 of Section 3. The 
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sequential gradient-restoration algorithm was employed in 
primal form (Refs. 16-18). This is a first-order algorithm 
which generates a sequence of feasible solutions, each characterized 
by a lower value of the performance index (10) . The numerical 
results are shown in Tables 1-2 and Fig. 2. 

Table 1 presents the following quantities: the entry values 
of the altitude, the velocity, and the path inclination; the 
exit values of the altitude, the velocity, and the path inclination; 
the minimum altitude, the peak heating rate, and the flight time; 
the characteristic velocity components and the total characteristic 
velocity . 

Table 2 compares the values of the total characteristic 
velocity for three trajectories: the optimal AOT trajectory, 
computed accounting for the control constraint (4); the ideal 
optimal AOT trajectory, computed disregarding the control 
constraint (4); this is the so-called grazing trajectory; and 
the optimal space trajectory; this is the so-called Hohmann 
transfer trajectory. 

IS 

Figure 2 contains four parts: the altitude h versus the 
time t (Fig. 2A) ; the velocity V versus the time t (Fig. 2B) ; 
the path inclination y versus the time t (Fig. 2C) ; and the lift 
coefficient C versus the time t (Fig. 2D) . 

J_i 

i l 

From Tables 1-2 and Fig. 2, the following comments arise: 

(i) The optimal trajectory includes two branches: a 
relatively short descending flight branch (branch 1) and a long 
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ascending flight branch (branch 2) . As the HEO radius increases, 
the minimum altitude of the optimal trajectory decreases, implying 
that a deeper penetration into the atmosphere is required to 
ensure the proper amount of velocity depletion. 

(ii) Velocity depletion takes place along the entire 
atmospheric trajectory, but is concentrated mostly in the terminal 
part of branch 1 and the beginning part of branch 2. 

(iii) The path inclination increases rapidly from the entry 
value (a few degrees negative) to zero value in branch 1; and 

it increases slowly from zero value to the exit value (a fraction 

of a degree positive) in branch 2. 

(iv) The lift coefficient profile is nearly independent 
of the HEO radius. In branch 1, the lift coefficient decreases 
rapidly from the upper bound value to nearly the lower bound value; 
in branch 2, the lift coefficient stays near the lower bound value. 

(v) For all the transfer cases studied, the characteristic 
velocity of the optimal trajectory is close to that of the ideal 
grazing trajectory. The excess expenditure in AV characterizing 
the optimal trajectory vis-a-vis the ideal grazing trajectory is 
1.5% for Case 1, is 0.9% for Case 2, and is 0.6% for Case 3 
(GEO-to-LEO transfer) . 

(vi) As the HEO radius increases, the advantage of the optimal 
trajectory vis-a-vis the Hohmann trajectory increases. The 
relative saving in AV is 52.2% in Case 1, is 57.8% in Case 2, and 


is 61.5% in Case 3. 
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4.3. Guidance Implications . Consider the altitude-path 
inclination domain and, with reference to Fig. 2, regard the 
altitude profile h = h(t) and the path inclination profile 
y = y(t) as parametric representation of the trajectory, the 
time t being the parameter. Upon elimination of the time, one 
obtains the path inclination-altitude relation y = y(h). Then, 
this relation can be rewritten in the normalized form 0 = 0 (n) , 
where 0 and n denote normalized variables defined as follows: 

6 = Y/|y 0 I , branch 1 , (11a) 

0 - Y/| Y Q | , branch 2, (lib) 

and 


n = (h Q - h)/(h Q - h min ) , branch 1, 

n = (h - h . )/(h n - h • ) / branch 2. 

1 nun'' 0 “mm 


(12a) 

(12b) 


The normalized path inclination-altitude relation 0 = 0(n) is 
plotted in Fig. 3, which contains two parts: the descending flight 
branch (Fig. 3A) and the ascending flight branch (Fig. 3B) . For 
branch 1, the normalized path inclination is nearly a linear function 
of the normalized altitude, and its slope is relatively steep; 
for branch 2, the normalized path inclination is also nearly a 
linear function of the normalized altitude, but its slope is 
relatively shallow. These observations are the basis of the 
gamma guidance law described in Section 5. This is a two-stage 
guidance law, designed as follows: for branch 1, the gamma 
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guidance is a linear path inclination guidance; for branch 2, 
the gamma guidance is a constant path inclination guidance. 
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5 . Gamma Guidance Trajectories 

In the previous sections, optimal trajectories for coplanar 
AOT flight were determined. They include three phases: the 
preatmospheric phase, characterized by the velocity impulse 
AV qo at HEO and the fact that the apogee of the HEO-to-entry 
transfer orbit occurs at r Q0 ; the atmospheric phase, characterized 
by properties (i)-(vi) of Section 4; and the postatmospheric 
phase, characterized by the fact that AV-^ = 0 at atmospheric exit, 
by the velocity impulse AV^ at LEO, and by the fact that the 
apogee of the exit-to-LEO transfer orbit occurs at r ll' 

In this section, we develop gamma guidance trajectories for 
coplanar AOT flight under two basic requirements: (a) the gamma 

guidance trajectory (GGT) should be close to the optimal trajectory 
(OT) ; and (b) the gamma guidance trajectory should be simple, easy 
to implement, and reliable. 

For the preatmospheric phase and the postatmospheric phase, 
we control the guidance trajectory via the velocity impulses 
AV q0 and AV^, albeit with slightly different values from those 
of the optimal trajectory, in order to increase reliability. In 
addition, at atmospheric exit, we introduce into the guidance 
trajectory the velocity impulse AV^ in order to compensate for 
previous velocity errors and ensure that the specified apogee 
can be reached. 

For the more complicated atmospheric phase, we structure 
the guidance trajectory so as to reproduce the two-branch 
geometry of the optimal trajectory: a relatively short descending 
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flight branch (branch 1) ; and a long ascending flight branch 
(branch 2) . To achieve the above properties, we employ a feedback 
control scheme, and we adjust the lift coefficient according 
to a two-stage gamma guidance law. For branch 1, the gamma 
guidance is a linear path inclination guidance; for branch 2, 
the gamma guidance is a constant path inclination guidance. We 
investigate the proper selection of four guidance parameters, 
namely: the initial path inclination y Q ; the target altitude h T 
of branch 1; the switch velocity Vg from branch 1 to branch 2; 
and the target path inclination y T of branch 2. 

5.1. Preatmospheric Phase . Initially, the spacecraft is in 
a high Earth orbit of radius r^. To deorbit, the following 
velocity impulse is applied: 


av 00 = /<r a /r 00 )v a 


' (r a /r 00 )V 0 COSY 0' 


(13a) 


with 

v 0 - v a /[2(r 00 2 - r 00 r a )/(r 00 2 - r \os 2 a ) \ . (13b) 

This enables the spacecraft to enter into an elliptical transfer 
orbit leading from HEO exit to atmospheric entry. In Eqs.(13), 
r rtn is the HEO radius, r is the radius of the outer edge of the 
atmosphere, and V is the circular velocity at r = r . Because 
r , are constant and r^^ is given, Eqs . (13) imply that 

v 0 = V Q (y 0 ) and AV Q0 = AV Q0 (y 0 ). Hence, the selection of the 
entry angle y^ determines uniquely both the entry velocity Vq 
and the initial velocity impulse AV QC) . 
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5.2. Atmospheric Phase . The atmospheric phase includes 

the descending flight branch (branch 1) and the ascending flight 
'branch (branch 2). For both branches, a gamma guidance scheme 
is implemented in feedback control form. The switch from 
branch 1 to branch 2 is regulated by the switch velocity , 
to be selected appropriately. 

Descending Flight Branch . For branch 1, the gamma guidance 
is a linear path inclination guidance, which is implemented in 
the following feedback control form: 


C L “ C L ( h / v / Y ) = -K x (Y - Y ) , 


Y = Y 0 (h - h T )/(h Q - h T ) , 


C_ £ C < C TK . 
La L — Lb 


(14a) 

(14b) 

(14c) 


Here, C is the instantaneous lift coefficient and C is the 
L L 

nominal lift coefficient; C T and C , are the lower and upper 
bounds for the lift coefficient; y is the instantaneous path 
inclination, y is the nominal path inclination, and Yq is the 
entry path inclination; h is the instantaneous altitude, h^ is 
the entry altitude, and h^ is the target altitude; and is 
the gain coefficient for path inclination error. With reference 
to (14) , the following remarks are pertinent. 

(a) The nominal lift coefficient C L is computed with 
Eq. (lc) under the assumption of near-equilibrium conditions. Upon 
setting y = 0, invoking Eqs. (2)-(3), and observing that p = p (h) , 
r = r(h), and g = g(h), we obtain the relation 
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(1/2 ) C^p (h) SV 2 + m[V 2 /r(h) - g(h)]cosy = 0, (15a) 

which admits the solution 

C L = 2m[g(h) - V 2 /r (h) ] cosy/p (h) SV 2 , (15b) 

which has the form 

C L = C L (h,V,y) . (15c) 

(b) The target altitude h^ should be in a proper range. 

If h,p is too high, the guidance trajectory is flat, resulting in 
an early exit from the atmosphere ; because not enough velocity 
is depleted in the atmospheric pass, this results in large values 
of the characteristic velocity. If h T is too low, the guidance 
trajectory is too steep, resulting in deeper penetration of 
the atmosphere; therefore, the peak heating rate is too high. 

In practice, the target altitude should be selected so as to 
be below, but close to, the minimum altitude of the OT. One 
option is to prescribe directly h ; another option is to prescribe 
indirectly h T , based on the selection of the target lift 
coefficient C^. The second option is preferable, because this 
procedure compensates automatically for variations of the atmospheric 
density with respect to the standard density. 

Next, observe that y = 0 at the target altitude. Then, 
with reference to Eq. (15b) , replace C L , h, V with C^ T , h T , V T , 
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where the subscript T denotes target value. This yields the 

relation 

C LT = 2m[g(h T ) - /r (h T ) ] /p (h T ) SV^, (16a) 

where 


V T = 0 . 25V q + 0.75V . (16b) 

Here, is the entry velocity of the OT and is the exit velocity 
of the OT. Because V , are known, V T is known. Hence, Eqs. (16) 
yield the one-to-one relation 


'LT 


C LT ' 


whose inverse form is 


(17a) 


h T =h T (C LT ) * (17b) 

Therefore, prescribing the target lift coefficient C is 

Li 1 

equivalent to prescribing the target altitude h,^. At any rate, 

the choice of C should be consistent with Ineq. (4) ; in 
-L 1 

addition, there should be enough margin on both the lower 
bound side and the upper bound side for control. 

(c) The gain coefficient for path inclination error is given 
by 


K x = p*/p, (18) 

where p = p (h) is the air density at the altitude h and p* = p (h*) 
is the air density at the reference altitude h* = H/3 = 40 km. 
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This particular form of the gain coefficient is justified by the 
need for a more energetic control response at higher altitudes 
and a gentler control response at lower altitudes. 

(d) Finally, it must be observed that the feedback form 
(14) of the gamma guidance has two important properties: (a) 

it is strongly stable at the target altitude; this is because 
y < 0 if h > h T# while y > 0 if h < h T ; and (b) it avoids 
overshooting and undershooting of the target altitude, since y 
varies smoothly between the entry value y = y Q and the target 
altitude value y = 0. 

Ascending Flight Branch . For branch 2, the gamma guidance 
is a constant path inclination guidance, which is implemented 
via the following feedback control form: 


■ C L (h,V,y) = -K 2 (y - y) , 

(19a) 

y t . 

(19b) 

± C L 1 C Lb- 

(19c) 


Here, y T denotes the target path inclination and is the gain 
coefficient for path inclination error. With reference to (19) , 
the following comments are pertinent. 

(a) The target path inclination y T should be in 
a proper range. If y T is too small, the exit from the atmosphere 
might become physically impossible. If y T is too large, the 
characteristic velocity component AV^ might become too large. 
Finally, the selection of y T should be consistent with the 
selection of the switch velocity V„ , so as to ensure that the 
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desired apogee can be reached after the atmospheric pass. 

(b) The gain coefficient for path inclination error is 
given by 


K 2 = p*/p, (20) 

where p = p (h) is the air density at the altitude h and p* = p(h*) 
is the air density at the reference altitude h* = H/3 = 40 km. 

Once more, this particular form of the gain coefficient is 

justified by the need for a more energetic control response 
at higher altitudes and a gentler control response at lower 
altitudes . 

5.3. Postatmospheric Phase . The postatmospheric phase 
includes two velocity impulses: a velocity impulse AV^ at 
atmospheric exit and a velocity impulse AV^ at LEO entry. 

Atmospheric Exit . The velocity impulse at atmospheric exit 
is determined with the relation 


AV 1 ^1 V 1 ' 


(21a) 


where 


V 1 = V [2(r U 2 - r ll r a )/(r ll 2 - r a 2 cos 2 Yl )]. (21b) 

This velocity impulse is essential for the GGT in order to 
compensate for previous velocity errors. In Eqs. (21), r^. is 
the LEO radius, r^ is the radius of the outer edge of the atmosphere, 
V a is the circular velocity at r = r Q , y^^ is the exit path 
inclination, is the exit velocity prior to the velocity impulse, 
and V 1 is the exit velocity after the application of the velocity 
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impulse. Note that r , are constant, is given, and 

arc measured in actual flight. 

LEO Entry. After the velocity impulse (21) is applied, the 


spacecraft enters into an elliptical transfer orbit leading from 
atmospheric exit to LEO entry. This elliptical transfer orbit 
is such that its apogee occurs at r^. At this point, the 
velocity impulse AV^ is applied so as to achieve circularization 
into LEO. Specifically, AV^ is determined with the relation 


AV 11 = /<r a /r ll )V a ' V ll- 


(22a) 


In Eq. (22a) , r is the radius of the outer edge of the atmosphere, 

V is the circular velocity at r = r , r . is the LEO radius, 
a an 

and is the velocity at LEO entry prior to the velocity impulse. 

Note that r v are constant, r. . is given, and V n , is measured 
a' v a 11 11 

in actual flight. However, in the computer simulations, the 
measurement of V.. is replaced by the relation 


V 11 = (r a /r ll )V l COSY l' 


(22b) 


which arises from angular momentum conservation applied to the 
exit-to-LEO transfer orbit. In Eq . (22b), V is supplied by 

Eq. (21b) . 

5.4. Guidance Parameters . As explained in the previous 
sections, the behavior of the guidance scheme described by 
Eqs . (13)- (22) depends on four parameters: the entry path 

inclination the target altitude h T ; 

and the target path inclination y . 


the switch velocity V g ; 
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(a) The entry path inclination y should be in a proper 
range. If |Yq| is too small, the guidance trajectory is too flat, 

( resulting in an early exit from the atmosphere; because not 
enough velocity is depleted in the atmospheric pass, this results 
in large values of the characteristic velocity. If |Yq| is 
too large, the guidance trajectory is too steep, resulting in 
deeper penetration of the atmosphere; therefore, the peak heating 
rate is too high. 

Good results are obtained by choosing |Yq| of the GGT to 
be somewhat larger than ly^of the OT. Note that the entry 
path inclination y^ determines automatically the entry velocity 
Vq and the HEO propulsive impulse AV^; see Eqs. (13). 

(b) The target altitude h T should not be too high nor 
too low. If h is too high, the guidance trajectory is flat, 
resulting in an early exit from the atmosphere; because not 
enough velocity is depleted in the atmospheric pass, this results 
in large values of the characteristic velocity. If h T is too low, 
the guidance trajectory is too steep, resulting in deeper 
penetration of the atmosphere; therefore, the peak heating rate 

is too high. 

Good results are obtained by choosing h T of the GGT to be 

somewhat lower than h . of the OT. As explained in Section 5.2, 

min 

this goal is achieved indirectly by prescribing the target lift 

coefficient C T corresponding to near-equilibrium level flight. 

J_i X 

(c) The switch velocity V should also be in a proper range. 
The best choice of Vg should be such that the exit velocity of 
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the GGT is close to the exit velocity of the OT. If this is 

done, the velocity impulse at atmospheric exit AV-^ is small. 

In practice, the switch velocity V of the GGT is chosen 

□ 

to be a weighted average of the entry velocity V Q of the OT and 
the exit velocity of the OT, 

v s = V i + A < v 0 " v i> ' (23a) 

where the dimensionless constant A is to be selected appropriately. 

Equation (23a) means that the switch velocity V is compatible with 

the target altitude h . Indeed, the higher the entry velocity 

V n , the lower the minimum altitude h . of the OT, hence 
u min 

the lower the target altitude h T of the GT; to guarantee 
the same value of the exit velocity , a higher switch velocity 
Vg is needed for the GT, which is the case with Eq. (23a) . 

(d) The target path inclination y T should be in a proper 
range. If y is too small, the exit from the atmosphere might 
become physically impossible. If y T is too large, the characteristic 
velocity component A\A^ might become too large. 

In practice, the target path inclination y T of the GGT is 
chosen to be some fraction of the exit path inclination y^ of 
the OT, 

Y T = B V ( 23b ) 


where the dimensionless constant B is to be selected appropriately. 
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5.5. Numerical Results . To investigate the effect of 
th° four guidance parameters y Qf h T , V , y t on the GGT, a 
large number of combinations were tried, and the corresponding 
guidance trajectories were computed. Typical combinations of 
guidance parameters are shown in Table 3 for Case 1, Case 2, 
and Case 3. These combinations were selected because the 
resulting GGT is geometrically close to the corresponding 
OT , so as to retain the good features of the OT concerning 
the total characteristic velocity 


AV 


AV 


00 


+ AV 1 + AV 


(24) 


and the peak heating rate. 

More detailed information about the GGT is shown in Table 4. 
This table presents the following quantities: the entry values 
of the altitude, the velocity, and the path inclination; the exit 
values of the altitude, the velocity, and the path inclination; 
the minimum altitude, the peak heating rate, and the flight time; 
the characteristic velocity components and the total characteristic 
velocity. 

Table 5 compares the GGT and the OT from the point of 
view of the total characteristic velocity AV. 

Tables 6-8 present a more detailed comparison of the GGT and 
the OT from the point of view of these quantities: the entry 
values of the altitude, the velocity, and the path inclination; 
the exit values of the altitude, the velocity, and the path 
inclination; the minimum altitude, the peak heating rate, and 
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the flight time; the characteristic velocity components and 
the total characteristic velocity. 

Finally, Figs. 4-6 present a graphical comparison of the GGT 
and the OT in terms of the altitude profile h(t), the velocity 
profile V ( t ) , the path inclination profile y(t), and the lift 
coefficient profile C (t) . 

From Tables 3-8 and Figs. 4-6, the following comments arise. 

(i) The altitude profile of the GGT is close to the 
altitude profile of the OT and retains the two-branch geometry 

of the OT: a relatively short descending flight branch (branch 1); 
and a long ascending flight branch (branch 2). As the HEO radius 
increases, the minimum altitude of the GGT decreases, consistently 
with the behavior of the OT. 

(ii) The velocity profile of the GGT is close to the 
velocity profile of the OT. Velocity depletion takes place 
along the entire atmospheric trajectory, but is concentrated 
mostly in the terminal part of branch 1 and the beginning part 
of branch 2 . 

(iii) The path inclination profile of the GGT is close 
to the path inclination profile of the OT. In branch 1, the 
path inclination increases almost linearly from the entry value 
to zero value; in branch 2, the path inclination increases slowly 
from zero value to the exit value. 

(iv) The lift coefficient profile of the GGT is close 
to the lift coefficient profile of the OT. In branch 1, the 
lift coefficient decreases rapidly from the upper bound value 

to the lower bound value; in branch 2, the lift coefficient stays 
near the lower bound value. 



30 


AAR-246 


(v) The characteristic velocity of the GGT is quite 

close to the characteristic velocity of the OT. The relative 
increase in AV is 1.3% for Case 1, 0.7% for Case 2, and 0.1% 
for Case 3. This small increase is almost exclusively due to 

the fact that AV^ / 0 for the GGT, while AV^ = 0 for the OT. 

In other words, while the OT exits from the atmosphere with the 
correct values of and y^, this is not the case with the GGT. 

Hence, at atmospheric exit, a nonzero, albeit small, velocity 
impulse is needed for the GGT. 

(vi) The peak heating rate (PHR) of the GGT is somewhat higher 
than the peak heating rate of the OT. The relative increase in 

PHR is 9.4% for Case 1, 4.7% for Case 2, and 7.7% for Case 3. 

The higher peak heating rate of the GGT vis-a-vis the OT is 

due to the larger value of |Yq| and the corresponding lower target 

altitude, which in turn cause a deeper penetration into the 
atmosphere. This increase in |y | and the corresponding decrease 
in h^ cannot be avoided, because of the need for stability of 
the gamma guidance trajectory. 
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6 . Modified Gamma Guidance Trajectories 

In real AOT flights, there are dispersion effects 
arising from navigation errors, variations of the atmospheric 
density, and uncertainties in the aerodynamic coefficients. 

Navigation errors refer to the space portion of AOT flights and 
induce errors in the entry path inclination; density variations 
are due to such factors as latitude, season, time of the day or 
the night, solar activity or are due to lack of sufficient knowledge 
of a particular planetary atmosphere (Mars) ; uncertainties in 
the aerodynamic coefficients arise because wind tunnel tests might 
not simulate precisely the combination of high speeds and low 
densities characterizing AOT flights or arise because CFD schemes^ 
might not account precisely for all of the physical factors involved. 

While the gamma guidance scheme of Section 5 yields a 
trajectory close to the optimal trajectory in the absence of 
parameter dispersion effects, this scheme is not sufficiently 
robust with respect to large parameter dispersion. For the 
sake of discussion , let unprimed quantities denote standard 
values; let primed quantities denote dispersed values; and 
let the following dispersion factors be defined: 


F 

YO 


= W 


(25a) 


F p = p' (h)/p(h) , 


(25b) 


CFD is an acronym for computational fluid dynamics. 
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F — C ' /C 

CDO DO 7 DO' 


K 


= K'/K, 


(25c) 
( 25d) 


F CLR = (C Lb - C ia> / (C Lb - C La>- 


(25e) 


Here, F is the entry path inclination factor; F is the density 
TO P 

factor; F^ d q is the zero-lift drag factor; F^ is the induced 

drag factor; and F is the lift range factor. If 

LLK 


there are no dispersion effects, 


F 

TO 


If 



F 

CDO 


If 


K 


= 1 , 


CLR 


= 1 . 


(26) 


However, if there are dispersion effects, one or more of the 
above factors might be different from unity. 

As an example, consider the entry path inclination factor 
(25a) in connection with Case 3; if F = 0.93, corresponding 
to AYq = +0.3, the GGT skips out of the atmosphere with little 
velocity dissipation; in turn, this results in a large value of 
the characteristic velocity. As another example, consider the 
density factor (25b) , and assume that its value is constant, 
independent of the altitude; if F^ = 1/3, the GGT skips out of 
the atmosphere. As a third example, consider the induced drag 
factor (25d) ; if F p = 1/2, the GGT skips out of the atmosphere; 
if Fp = 2, the GGT penetrates too deeply into the atmosphere, 
ultimately leading to a crash. 

There are two ways for improving the stability of the gamma 
guidance scheme with respect to dispersion effects: 
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(i) to lower the target altitude; (ii) to adjust the switch 
velocity and the target path inclination by means of a predictor- 
corrector algorithm. The resulting trajectory is called modified 
gamma guidance trajectory (MGGT) . 

6.1. Target Altitude and Entry Path Inclination . From the 
point of view of the total characteristic velocity and the peak 
heating rate, the target altitude should be as high as possible; 
hence, it should be as close as possible to the minimum altitude 
of the optimal trajectory. While the GGT of Section 5 achieves 
this goal, it is marginally stable with respect to dispersion 
effects . 

With the MGGT, the goal is to increase the stability with 
respect to dispersion effects, while keeping the total characteristic 
velocity and the peak heating rate within reasonable range. The 
increased stability is achieved by lowering the target altitude 
h , hence increasing the target lift coefficient C . To achieve 
compatibility between the target altitude and the entry path 
inclination y , steeper values of the entry path inclination might 
be desirable. 

6.2. Switch Velocity and Target Path Inclination . Because 
of parameter dispersion effects, the ideal condition AV^ = 0 
cannot be achieved . However , AV, can be kept small by adjusting 
the switch velocity and the target path inclination in such a way 
that AV^ is less than some threshold value, for instance, 

AV^ < 0.03 km/sec. This can be achieved via the predictor-corrector 
algorithm described below. 
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6.3. Predictor-Corrector Algorithm . This algorithm includes 

five steps. The description of these steps is given below 

with particular reference to Case 3, GEO— to-LEO transfer. 

1 

Step 1 . Assume that the spacecraft is at the end of branch 1 
with switch altitude and switch path inclination given by 


h 


S 




0 . 


Assume that the switch velocity is 


(27) 


v s = v i + A(V Q - V 1 ) (28a) 

and that the target path inclination of branch 2 is 

Y T = BY r < 28b > 

Here, is the entry velocity of the optimal trajectory; is 
the exit velocity of the optimal trajectory; is the exit path 
inclination of the optimal trajectory; and A, B are dimensionless 
constants to be specified appropriately. 

The predictor-corrector algorithm starts by assuming that 

V = 8.40 km/sec, y t = 0.15 deg, (29a) 

corresponding to 

A = 0.25, B = 0.50. (29b) 

Step 2. For branch 2, integrate Eqs. (l)-(3) in forward 
time using the feedback control form (19) of the constant path 
inclination guidance. Determine the exit values of the 

guidance trajectory. 
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Step 3 . Verify whether the inequality 

AV^ <_ 0.03 km/sec (30) 

is satisfied. If this is the case, accept the assumed pair 
(A,B) , hence the assumed pair (Vg, Y T ) ; the predicted-corrector 
procedure is terminated. Otherwise, go to either Step 4 (if V-^ 
must be increased) or Step 5 (if must be decreased). 

Step 4 . This step must be performed if compliance with 
Ineq. (30) requires that be increased. It consists of two 
successive one-dimensional searches, first in the A-direction, 
while keeping B = const (Step 4A) , and then in the B-direction, 

while keeping A = const (Step 4B) . 

Step 4A . For B = const, hence y T = const, subject the 
parameter A to increments AA corresponding to increments 
AVg = 0.10 km/sec. For each value of A, hence for each value of 

Vg, execute Step 2. The search in the A-direction stops if 

Ineq. (30) is satisfied simultaneously with V < 9.10 km/sec; in 
such a case, the assumed (A,B) pair is accepted, and the predictor- 
corrector procedure is terminated. However, if the value 
Vg = 9.10 km/sec is reached without satisfaction of Ineq. (30), 
then go to Step 4B. 

Step 4B. For A = const, corresponding to V„ = 9.10 km/sec, 
subject the parameter B to increments AB corresponding to 
increments Ay^ = 0.01 deg. For each value of B, hence for each 
value of y t , execute Step 2. The search in the B-direction stops 
if Ineq. (30) is satisfied; in such a case, the assumed (A,B) 
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pair is accepted, and the predictor-corrector procedure is 
terminated . 

Step 5 . This step must be performed if compliance with 
Ineq. (30) requires that V 1 be decreased. It consists of a 
single one-dimensional search in the A-direction, while keeping 
B = const. 

For B = const, hence y = const, subject the parameter A 
to decrements AA corresponding to decrements AV g = -0.05 km/sec. 

For each value of A, hence for each value of V , execute Step 2. 

D 

The search in the A-direction stops if Ineq. (30) is satisfied; 
in such a case, the assumed (A,B) pair is accepted, and the 
predictor-corrector procedure is terminated. 

6.4. Numerical Results . Here, we present a comparative 
study of the GGT and the MGGT for Case 3, GEO-to-LEO transfer. 

This study is done in the absence of dispersion effects, that is, 
under assumptions (26) . Dispersion effects are analyzed 
systematically in Section 7. 

For both the GGT and the MGGT, Table 9 presents typical 
combinations of the guidance parameters. 

For both the GGT and the MGGT, Table 10 shows the following 
quantities: the entry values of the altitude, the velocity and 
the path inclination; the exit values of the altitude, the velocity, 
and the path inclination; the minimum altitude, the peak heating 
rate, and the flight time; the characteristic velocity components 
and the total characteristic velocity. 

Finally, Fig. 7 presents a graphical comparison of the GGT 
and the MGGT in terms of the altitude profile h(t), the velocity 



37 


AAR-246 


profile V(t), the path inclination profile y(t), and the lift 
coefficient profile C (t) . 

Xj 

From Tables 9-10 and Fig. 7, following comments arise. 

(i) The altitude profile of the MGGT lies below the altitude 
profile of the GGT. 

(ii) The velocity profile of the MGGT is relatively close 
to the velocity profile of the GGT. 

(iii) The path inclination profile of the MGGT is relatively 
close to the path inclination profile of the GGT. 

(iv) The lift coefficient profile of the MGGT is quite 
different from the lift coefficient profile of the GGT in that 
an ampler lift coefficient margin now exists with respect to the 
lower bound value. 

(v) The characteristic velocity of the MGGT is quite 
close to the characteristic velocity of the GGT. The relative 
increase in AV is 0.4%. 

(vi) The peak heating rate of the MGGT is somewhat higher 
than the peak heating rate of the GGT. The relative increase 

is 23. 9%. The higher peak heating rate of the MGGT is due to 
deeper penetration into the atmosphere (by 3.1 km), which in turn 
is tied to the larger value of |Yq| and the lower value of h T . 

This increase in |y | and decrease in h^ constitute the price 
which must be paid for the increased stability of the modified 
gamma guidance trajectory. 
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7 . Parameter Dispersion Effects 

As pointed out in Section 6, parameter dispersion eirects 
are important in real AOT flights; therefore, they must be 
considered in the evaluation of the merits of a particular 
guidance scheme. In this section, we consider the modified 
gamma guidance scheme of Section 6 and explore the behavior of 
the MGGT vis-a-vis deviations from unity of the dispersion 
factors (25) . 

With reference to Case 3, GEC-to-LEO transfer, the results 
are shown in Tables 11-15 and Fig. 8. They demonstrate that the 
MGGT is capable of executing safely the atmospheric pass in the 
following range of values of the dispersion factors: 


0.91 

< 

F 

< 

1.44 , 

(31a) 



Y0 




0.33 

< 

F 

P 

< 

10.00, 

(31b) 

0.10 

< 

F 

< 

5.00, 

(31c) 



CD0 




0.10 

< 

F^ 

< 

5.00, 

Old) 



K 




0.33 

< 

F„ T ^ 

< 

3.00. 

(31e) 


— 

CLR 

— 




A detailed discussion follows. 

7.1. Entry Path Inclination Factor . See Table 11 and 
Fig. 8A. The MGGT is stable for values of the entry path 
inclination factor in the range 

YO 


0.91 < F 


< 1.44, 


(32a) 
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corresponding to 

-4.1 >_ Yq -6.5 deg, (32b) 

+ 0.4 >_ AYq >_ -2.0 deg. (32c) 

In the above range, the target altitude, the switch velocity, 

and the target path inclination are constant. It must be noted 

that, if = 0.89, corresponding to Yq = - 4.0 deg and AYq = +0.5 deg, 

the MGGT skips out of the atmosphere. On the other hand, if 

F = 1.44, corresponding to Yq = - 6.5 deg and AYq = -2.0 deg, 

the MGGT undershoots the target altitude. This situation is due 

to the fact that only a simple proportional feedback control is 

used in this paper and can be alleviated by employing more 

sophisticated forms of feedback control. 

7.2. Density Factor . See Table 12 and Fig. 8B. The MGGT 
is stable for values of the density factor in the range 

0.33 < F < 10.00. (33) 

- p - 

In the above range, the entry path inclination, the switch 

velocity, and the target path inclination are constant. Both 

the target altitude and the peak heating rate increase as F^ 

increases. It is clear from (33) that the modified gamma 

guidance can tolerate large density increases better than large 

density decreases. In particular, if F^ = 0.25, the MGGT skips 

out of the atmosphere with a large increase in characteristic velocity. 
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7.3. Zero-Lift Drag Factor . See Table 13 and Fig. 8C. 

The MGGT is stable for values of the zero-lift drag factor in 
the range 

0.10 < F cd0 < 5.00. (34) 

In the above range, the entry path inclination and the target 

altitude are constant. As F increases, the switch velocity 

CD0 

increases, while the flight time decreases considerably. The 
increase in switch velocity is due to the fact that, as 
increases, more energy is depleted for ascending from 
the target altitude to the atmospheric exit altitude. It is 
clear from (34) that the modified gamma guidance is able to 
tolerate both large increases and large decreases in zero-lift drag. 

7.4. Induced Drag Factor . See Table 14 and Fig. 8D. The 
MGGT is stable for values of the induced drag factor in the range 

0.10 < F < 5.00. (35) 

— K — 

In the above range, the entry path inclination and the 

target altitude are constant. As F^ increases/ 

the switch velocity increases, while the flight time decreases 

considerably. The increase in switch velocity is due to the 

fact that, as F K increases, more energy is depleted for 

ascending from the target altitude to the atmospheric exit altitude. 
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It is clear from (35) that the modified gamma guidance is able 
to tolerate both large increases and large decreases in induced drag. 

7.5. Lift Range Factor . See Table 15 and Fig. 8E. The 
MGGT is stable for values of the lift range factor such that 

0.33 < F^ t _ < 3.00. (36) 

In the above range, the entry path inclination and the target 
path inclination are constant. As F CLR increases, the target 
altitude increases, the switch velocity decreases, and the target 
lift coefficient, which is negative, becomes larger in modulus; 
indeed, C is proportional to F . As F increases, the peak 
heating rate decreases considerably, and the flight time changes 
somewhat. Generally speaking, the MGGT is more able to tolerate 
large increases in F _ T than large decreases. In particular, 

C-LiK 

if F =0.25, the MGGT skips out of the atmosphere with a large 
CLR 

increase in characteristic velocity. 
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8 . Conclusions 

This paper is concerned with the optimization and 
guidance of trajectories for coplanar, aeroassisted orbital 
transfer from high Earth orbit to low Earth orbit. The 
following major conclusions are obtained: 

(i) The optimal trajectories include two branches: a 
relatively short descending flight branch (branch 1) and a long 
ascending flight branch (branch 2). In branch 1, the path 
inclination ranges from a few degrees negative to zero and is 
nearly a linear function of the altitude; in branch 2, the path 
inclination ranges from zero to a fraction of a degree positive 
and is a slowly varying function of the altitude. 

(ii) Gamma guidance trajectories are developed so as to 
approximate the optimal trajectories in real time, while 
retaining the essential characteristics of simplicity, ease of 
implementation, and reliability. For the atmospheric pass, a 
feedback control scheme is employed and the lift coefficient 
is adjusted according to a two-stage gamma guidance law. For 
branch 1, the gamma guidance is a linear path inclination guidance; 
for branch 2, the gamma guidance is a constant path inclination 
guidance. By proper selection of four guidance parameters (the 
entry path inclination, the target altitude of branch 1, the 
switch velocity, and the target path inclination of branch 2), 

the gamma guidance trajectory can be made close to the optimal 
trajectory . 



43 


AAR-246 


(iii) Improvements in stability are possible via a modified 
gamma guidance, which differs from the gamma guidance in that 
the switch velocity and the target path inclination are adjusted 
in flight with a predictor-corrector algorithm; also, the target 
altitude is lower and the entry path inclination is steeper. 
Computer simulations show that the modified gamma guidance 
trajectory is quite stable with respect to dispersion 

effects arising from navigation errors, variations of the 
atmospheric density, and uncertainties in the aerodynamic 
coefficients . 

(iv) A byproduct of the parameter dispersion studies 
is the following design concept. For coplanar, aeroassisted 
or bital transfer, the lif t-range-to-weight ratio appears to play 

a more important role than the lift-to-drag ratio. This is because 
the lif t-range-to-weight ratio controls mainly the minimum altitude 
(hence, the peak heating rate) of the guidance trajectory; on 
the other hand, the lift-to-drag ratio controls mainly the 
duration of the atmospheric pass of the guidance trajectory. 
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Table 1. Results 

for the optimal 

trajectories. 


Quantity 

Case 1 

Case 2 

Case 3 

Units 

h o 

120.0 

120.0 

120.0 

km 

v o 

9.040 

9.905 

10.310 

km/sec 

Y 0 

-3.034 

-3.893 

-4.204 

deg 

h l 

120.0 

120.0 

120.0 

km 

V 1 

7.844 

7.844 

7.844 

km/sec 

*1 

0.319 

0.319 

0.319 

deg 

^min 

79.50 

76.35 

75.36 

km 

PHR 

35.90 

59.61 

72.70 

W/crrt^ 

T 

2.147 

2.297 

2.347 

ksec 

o 

o 

> 

<a 

1.025 

1.445 

1.490 

km/ sec 

av i 

0.000 

0 .000 

0.000 

km/sec 

Av n 

0.025 

0.025 

0.025 

km/sec 

AV 

1.049 

1.470 

1.515 

km/ sec 
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Table 2. Comparison of characteristic velocities AV(km/sec). 


Trajectory 

Case 1 

Case 2 

Case 3 

Optimal 

1.049 

1.470 

1.515 

Grazing 

1.034 

1.457 

1.504 

Hohmann 

2.195 

3.486 

3.940 
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Table 

3. Parameters 

of gamma 

guidance trajectories. 

Quantity 

Case 1 

Case 2 

Case 3 

Units 

Y 0 

-3.300 

-4.100 

-4.400 

deg 

N 

78.25 

75.75 

74.32 

km 

v s 

8.700 

9.180 

9.400 

km/sec 

y t 

0.115 

0.113 

0.114 

deg 

c 

^LT 

-0.270 

-0.270 

-0.288 

- 

A 

0.716 

0.648 

0.632 



B 


0.361 


0.354 


0.357 
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Table 

4 . Results 

for gamma 

guidance trajectories 

• 

Quantity 

Case 1 

Case 2 

Case 3 

Units 

h o 

120.0 

120.0 

120.0 

km 

V 0 

9.039 

9.905 

10.310 

km/sec 

Y o 

-3.300 

-4.100 

-4.400 

deg 

h l 

120.0 

120.0 

120.0 

km 

V 1 

7.832 

7.834 

7.841 

km/sec 

*1 

0.187 

0.217 

0.300 

deg 

h 

78.29 

75.73 

74.28 

km 

mm 





PHR 

39 .29 

62.44 

78.27 

2 

W/cm 

T 

2.456 

2.527 

2.405 

ksec 

O 

o 

> 

<1 

1.026 

1.446 

1.490 

km/sec 

o 

< 

0.016 

0.013 

0.003 

km/sec 

Av n 

0.020 

0.021 

0.024 

km/ sec 


km/ sec 


AV 


1.063 


1.480 


1.517 
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Table 5. Comparison 

of characteristic 

velocities 

AV (km/sec ) . 


Trajectory 

Case 1 

Case 2 

Case 3 


OT 

1.049 

1.470 

1.515 


GGT 

1.063 

1.480 

1.517 
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Table 6 . 

Comparison of 

results, Case 1. 


Quantity 

OT 

GGT 

Units 

h 0 

120.0 

120.0 

km 

V 0 

9.040 

9.039 

km/ sec 

^0 

-3.034 

-3.300 

deg 

h l 

120.0 

120.0 

km 

V 1 

7.844 

7.832 

km/sec 

Y 1 

0.319 

0.187 

deg 

h . 

79.50 

78.29 

km 

min 




PHR 

35.90 

39.29 

2 

W/cm 

T 

2.147 

2.456 

ksec 

o 

o 

> 

<3 

1.025 

1.026 

km/sec 

Av i 

0.000 

0.016 

km/sec 

Av n 

0.025 

0.020 

km/sec 

AV 

1.049 

1.063 

km/sec 
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Table 7 . 

Comparison 

of results. Case 2. 


Quantity 

OT 

GGT 

Units 

h 0 

120.0 

120.0 

km 

v o 

9.905 

9.905 

km/sec 

Y 0 

-3.893 

-4.100 

deg 

h l 

120.0 

120.0 

km 

V 1 

7.844 

7.834 

km/sec 

Y 1 

0.319 

0.217 

deg 

^min 

76.35 

75.73 

km 

PHR 

59.61 

62.44 

W/cm^ 

T 

2.297 

2.527 

ksec 

t> 

< 

o 

o 

1.445 

1.446 

km/sec 

AV 1 

0.000 

0.013 

km/ sec 

Av ll 

0.025 

0.021 

km/sec 

AV 

1.470 

1.480 

km/sec 
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Table 8 . 

Comparison of 

results. Case 3. 


Quantity 

OT 

GGT 

Units 

h o 

120.0 

120.0 

km 

v o 

10.310 

10.310 

km/ sec 

Y o 

-4.204 

-4.400 

deg 

h l 

120.0 

120.0 

km 

V 1 

7.844 

7.841 

km/sec 

Y 1 

0.319 

0.300 

deg 

h . 

75.36 

74.28 

km 

min 




PHR 

72.70 

78.27 

W/cm^ 

T 

2.347 

2.405 

ksec 

o 

o 

> 

< 

1.490 

1.490 

km/sec 

1 — 1 
> 
< 

0.000 

0.003 

km/sec 

Av ll 

0.025 

0.024 

km/sec 

AV 

1.515 

1.517 

km/sec 
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Table 9. 

Parameters 

of modified gamma guidance 

trajectories, Case 3 

Quantity 

GGT 

MGGT 

Units 

Y 0 

-4.400 

-4.500 

deg 

h T 

74.32 

71.14 

km 

v s 

9.400 

8.400 

km/ sec 

y t 

0.114 

0.150 

deg 

C LT 

-0.288 

-0.180 

- 

A 

0.632 

0.225 

- 

B 

0.357 

0.470 

- 

F = 1, 

YO 

f p ■ 

F =1 F=1.F =" 

CDO ' K ' CLR 

L. 


1 
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Table 10. 

Results 

for modified 

gamma guidance trajectories. Case 3 

Quantity 

GGT 

MGGT 

Units 

h 0 

120.0 

120.0 

km 

v o 

10.310 

10.310 

km/ sec 

Y o 

-4.400 

-4.500 

deg 

h l 

120.0 

120.0 

km 

V 1 

7.841 

7.835 

km/sec 

Y 1 

0.300 

0.245 

deg 

^min 

74.28 

71.14 

km 

PHR 

78.27 

96 . 97 

W/cm^ 

T 

2.405 

3.012 

ksec 

Av oo 

1.490 

1.491 

km/ sec 

i— 1 
> 
< 

0.003 

0.011 

km/ sec 

Av n 

0.024 

0.022 

km/sec 

AV 

1.517 

1.523 

km/sec 

F =1, 

Y0 

F p = 1, 

F =1 

CDO 

f k - f clr = !• 
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iaoie ha . 


Effect of path inclination factor F 
modified gamma guidance trajectories 0 ’ 


Case 3, 


F 

Y0 

Y 0 

h m 


deg 

km 

0.91 

-4.1 

71.14 

0.96 

-4.3 

71.14 

1.00 

-4.5 

71.14 

1.11 

-5.0 

71.14 

1.22 

-5.5 

71.14 

1.44 

-6.5 

71.14 


deg 


0.225 0.470 
0.225 0.470 
0.225 0.470 

°* 225 0.470 
0.225 0.470 
0.225 0.470 


'S 

km/sec 


Yn 


'LT 


B 


8.4 


u.lO 

00 

* 

O 

1 

0.15 

-0.18 

0.15 

-0.18 

0.15 

-0.18 

0.15 

-0.18 

0.15 

-0.18 
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Table 11B. Effect of path inclination factor F , Case 3, 

YO 

modified gamma guidance trajectories. 


— 

f yo 

h o 

km 

v o 

km/ sec 

Y 0 

deg 

h l 

km 

V 1 

km/sec 

Y 1 

deg 

— 

0.91 

120.0 

10.310 

-4.1 

120.0 

7.835 

0.244 

— 

0.96 

120.0 

10.310 

-4.3 

120.0 

7.836 

0.245 


1.00 

120.0 

10.310 

-4.5 

120.0 

7.835 

0.245 


1.11 

120.0 

10.310 

-5.0 

120.0 

7.835 

0.244 


1.22 

120.0 

10.309 

-5.5 

120.0 

7.836 

0.245 


1.44 

120.0 

10.309 

-6.5 

120.0 

7.835 

0.243 
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Table lie. Effect of path inclination factor F , Case 3, 

YO ' 

modified gamma guidance trajectories. 


F 

YO 

h . 
mm 

km 

PHR 

W/cm^ 

T 

ksec 

AV nn 

00 

km/ sec 

AV 1 

km/sec 

AV 11 
km/ sec 

AV 

km/sec 

0.91 

71.14 

83.79 

2.986 

1.490 

0.011 

0.022 

1.523 

0.96 

71.14 

93.33 

3.010 

1.490 

0.011 

0.022 

1.523 

1.00 

71.14 

96.97 

3.012 

1.491 

0.011 

0.022 

1.523 

1.11 

69.67 

113.45 

2.982 

1.492 

0.011 

0.022 

1.525 

1.22 

66.43 

136.61 

2.927 

1.493 

0.011 

0.022 

1.526 

1.44 

60.91 

180.80 

2.810 

1.496 

0.011 

0.022 

1.529 
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Table 12A. Effect of density factor F , Case 3, 
modified gamma guidance trajectories. 
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Table 12B. Effect of density factor F , Case 3, 
modified gamma guidance trajectories. 


— 

F 

P 

h o 

km 

v o 

km/sec 

Y 0 

deg 

h l 

km 

V 1 

km/sec 

Y 1 

deg 


0.33 

120.0 

10.310 

-4.5 

120.0 

7.824 

0.116 


0.50 

120.0 

10.310 

-4.5 

120.0 

7 .828 

0.175 


1.00 

120.0 

10.310 

-4.5 

120.0 

7.835 

0.245 


2.00 

120.0 

10.310 

-4.5 

120.0 

7.842 

0.273 


3.00 

120.0 

10.310 

-4.5 

120.0 

7.845 

0.276 

— 

5.00 

120.0 

10.310 

-4.5 

120.0 

7.849 

0.267 


10.00 

120.0 

10 .310 

-4.5 

120.0 

7.853 

0.241 
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Table 12C. Effect of density factor F^, Case 3, 
modified gamma guidance trajectories. 


F 

P 

h . 
min 

km 

PHR 
W/ cra^ 

T 

ksec 

Av oo 

km/ sec 

AV 1 

km/ sec 

iv n 

km/sec 

AV 

km/ sec 

0.33 

62.98 

87.66 

3.670 

1.491 

0.025 

0.019 

1.535 

~ 0.50 

66.13 

91.36 

3.366 

1.491 

0.020 

0.020 

1.531 

1.00 

71.14 

96.97 

3.012 

1.491 

0.011 

0.022 

1.523 

2.00 

75.52 

103.99 

2.773 

1.491 

0.003 

0.023 

1.517 

_ 3.00 

77.67 

109.09 

2.660 

1.491 

0.000 

0.023 

1.514 

5.00 

80.08 

116.50 

2.535 

1.491 

0.003 

0.023 

1.516 

"10.00 

83.11 

127.67 

2.376 

1.491 

0.007 

0.022 

1.519 
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Table 13A. Effect of zero-lift drag factor F^q, Case 3, 
modified gamma guidance trajectories. 



F 

CDO 

Y o 

h T 

v s 

y t 

C LT 

A 

B 

— 


deg 

km 

km/sec 

deg 





71.14 

8.20 

0.15 

00 
1 — l 

• 

0 

1 

0.144 

0.470 

i — I 
• 

rH 

r** 

8.25 

0.15 

-0.18 

0.165 

0.470 

71.14 

8.35 

0.15 

-0.18 

0.205 

0.470 

71.14 

8.40 

0.15 

1 

o 

• 

H* 

CD 

0.225 

0.470 

71.14 

8.50 

0.15 

00 
i — 1 

• 

0 

1 

0.266 

0.470 

71.14 

9.10 

0.21 

-0.18 

0.509 

0.658 
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Table 13B. Effect of zero-lift drag factor F^q, Case 3, 
modified gamma guidance trajectories. 

f cdo h o v 0 Y o h i V 1 Y 1 

km km/sec deg km km/sec deg 

0.10 120.0 10.310 -4.5 120.0 7.861 0.462 

0.20 120.0 10.310 -4.5 120.0 7.863 0.476 

0.50 120.0 10.310 -4.5 120.0 7.859 0.447 

1.00 120.0 10.310 -4.5 120.0 7.835 0.245 

2.00 120.0 10.310 -4.5 120.0 7.823 0.114 

5.00 120.0 10.310 -4.5 120.0 7.825 0.197 
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Table 13C. Effect of zero-lift drag factor f CD q' Case 3, 
modified gamma guidance trajectories. 



F 

CDO 

h . 
mm 

km 

PHR 

W/cm^ 

T 

ksec 

Av oo 

km/sec 

AV 1 

km/sec 

AV 11 

km/sec 

AV 

km/sec 


0.10 

71.14 

98.96 

4.268 

1.491 

0.025 

0.032 

1.547 

— 

0.20 

71.14 

98.73 

3.686 

1.491 

0.028 

0.033 

1.551 


0.50 

71.14 

98.05 

3.043 

1.491 

0.022 

0.031 

1.543 


1.00 

71.14 

96.97 

3.012 

1.491 

0.011 

0.022 

1.523 


2.00 

71.14 

94.98 

3.024 

1.491 

0.026 

0.019 

1.536 


5.00 

71.14 

90.00 

1.953 

1.491 

0.022 

0.020 

1.534 

— 
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Table 14A. Effect of induced drag factor F , Case 3, 

K 

modified gamma guidance trajectories. 


Y0 

deg 

hip 

km 

V S 

km/sec 

y t 

deg 

C LT 

A 

B 

■4.5 

71.14 

8.15 

0.15 

-0.18 

0.124 

0.470 

•4.5 

71.14 

8.15 

0.15 

-0.18 

0.124 

0.470 

•4.5 

71.14 

8.25 

0.15 

-0.18 

0.165 

0.470 

•4.5 

71.14 

8.40 

0.15 

-0.18 

0.225 

0.470 

■4.5 

71.14 

8.40 

0.15 

-0.18 

0.225 

0.470 

•4.5 

71.14 

9.10 

0.25 

-0.18 

0.509 

0.784 
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Table 14B. Effect of induced drag factor F , Case 3, 

K 

modified gamma guidance trajectories. 
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— 


Table 14C. 

Effect of 
modified 

induced drag factor F , Case 

IN. 

gamma guidance trajectories. 

3, 


— 

F 

K 

h . 
mm 

km 

PHR 

W/cm^ 

T 

ksec 

Av oo 

km/sec 

AV 1 

km/ sec 

iv n 

km/sec 

AV 

km/sec 

— 

0.10 

71.14 

103.63 

4.281 

1.491 

0.011 

0.021 

1.523 


0.20 

71.14 

102.66 

4.183 

1.491 

0.016 

0.020 

1.527 


0.50 

71.14 

100.19 

3.186 

1.491 

0.026 

0.031 

1.548 

— 

1.00 

71.14 

96.97 

3.012 

1.491 

0.011 

0.022 

1.523 


2.00 

71.14 

92.33 

3.072 

1.491 

0.027 

0.019 

1.537 


5.00 

71.21 

83.95 

1.646 

1.491 

0.029 

0.021 

1.541 
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Table 15A, 

. Effect of 
modified 

lift range factor F , Case 

LLK 

gamma guidance trajectories. 

3 , 



f clr 

Y 0 

h 

T 

V S 

y t 

C LT 

A 

B 

— 


deg 

km 

km/ sec 

deg 




— 

0.33 

-4.5 

62.98 

9.0 

0.15 

-0.059 

0.469 

0.470 


0.50 

-4.5 

66.13 

8.5 

0.15 

-0.090 

0.266 

0.470 


1.00 

-4.5 

71.14 

8.4 

0.15 

-0.180 

0.225 

0.470 

— 

2.00 

-4.5 

75.79 

8.4 

0.15 

-0.360 

0.225 

0.470 


3.00 

-4.5 

78.41 

8.4 

0.15 

-0.540 

0.225 

0.470 
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Table 15B. Effect of lift range factor F_ T _, Case 3, 

CLR 

modified gamma guidance trajectories. 
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Table 15C. Effect of lift range factor Case 3, 

modified gamma guidance trajectories. 



F 

CLR 

h . 
min 

km 

PHR 

W/cm^ 

T 

ksec 

Av oo 

km/sec 

AV 1 

km/sec 

Av n 

km/ sec 

AV 

km/sec 

— 

0.33 

61.50 

174.29 

2.390 

1.491 

0.029 

0.036 

1.555 

— 

0.50 

64.41 

157.37 

3.567 

1.491 

0.029 

0.018 

1.538 


1.00 

71.14 

96.97 

3.012 

1.491 

0.011 

0.022 

1.523 


2.00 

75.79 

62.44 

2.879 

1.491 

0.022 

0.019 

1.532 


3.00 

78.41 

49 . 08 

2.571 

1.491 

0.023 

0.019 

1.533 
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150 



FIG. 5A . TRAJECTORY COMPARISON , 
CASE2 . 


14 



FIG. 5B . TRAJECTORY COMPARISON. 
CASE2 . 
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FIG. 5C . TRAJECTORY COMPARISON, 
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FIG. 5D . TRAJECTORY COMPARISON, 
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O FGAMMR0=0.91 
A FGRMMR0=1.00 
+ FGRMMA0=1.44 


8R. MODIFIED GUIDANCE TRAJECTORIES. 
ENTRY GAMMA CHANGE. CASE3 . 


O FRH0=0 . 33 
* FRH0=1 .00 
+ FRH0=3 .00 


8B. MODIFIED GUIDANCE TRAJECTORIES. 
DENSITY CHANGE. CASE3 . 
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H(KM) 


O FCD0=0 . 1 
A FCDOrl.O 
+ FCD0=5.0 


FIG. 8C . MODIFIED GUIDANCE TRAJECTORIES , 
ZERO-LIFT DRAG CHANGE. CASE3 . 


H( KM ) 


O FK=0 . 1 
A FKrl.O 
+ FK=5.0 


FIG. 8D • MODIFIED GUIDANCE TRAJECTORIES. 
INDUCED DRAG CHANGE. CASE3 . 








